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ABSTRACT 


In order to establish a model which would more accurately predict 
discharge characteristics of a cell at any arbitrary current, and which 
would better fit the family of discharge curves, testing was done based 
on the Shepherd equation using a transparent lead-acid battery. 


A set of discharge curves at .6, 1.5, 3.6, and 5.4A, preceded by a 
3.6A preconditioning discharge to remove the "memory effect," was obtained. 
The battery was cycled continuously through deep discharges in an attempt 
to determine secular changes during the life of the battery. 


Several modifications to the Shepherd equation were proposed in 
order to obtain a better fit to discharge data. These included modeling 
the capacity as a function of current as Q=CI 1-n), Peukert's equation; 
modeling the ohmic polarization as a linear function of accumulated charge; 
and modeling the diffusion polarization as a function of accumulated charge 
alone and not current. 


Conclusions were drawn on the suitability of Shepherd's equation. 
The inclusion of the changed ohmic polarization term gave a good fit for 
individual curves which also suggested the change in the diffusion polari- 
zation term. Since a small current density approximation was used in 
Shepherd's equation, developing a model based on the original terms was 
recommended. Also recommended were the modification of a model to predict 
characteristics during a cycle with various discharge currents and further 
life cycle testing to determine the secular changes in the battery. 
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CHAPTER l 
INTRODUCTION 
1.1 BACKGROUND 

Predicting the transient and long-term characteristics 
of a cell under arbitrary charging and discharging histories 
is a problem confronting a system designer which could be 
overcome if there were an electrochemical model of the 
lead-acid cell. With this model, it would theoretically 
become possible to determine the optimum scheme for getting 
energy into and out of a battery. 

To what degree is the efficiency of a battery affected 
by different types of discharge? Should the designer select 
a pulse-type discharge over a steady discharge? What effect 
does this have on the system as a whole? To what degree does 
the age of the battery affect the efficiency? Low current or 
high? What exactly is the effect of temperature and current 
on the charging and discharging history of the arbitrary 
battery? Will current models handle the problem of the age 
of the battery? Does the model provide a fit merely to a 
Single discharge curve or for the family of curves? Does it 
account for the condition of the battery? What effect has 
gassing on the capacity for full charge or discharge? 
Questions such as these pose a great problem for the designer 
of a model. 


Naturally, there arenumerous other criteria which might 
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affect the optimum charge/discharge history of a battery, so 
it becomes very important to select a battery type to test 
which possesses characteristics which are relatively similar 
ferocher cells. 

The lead-acid cell was chosen for this investigation for 
several reasons. lLead-acid batteries are the most widely 
used batteries in the world, and for that reason, they 
represent the greatest commercial interest. Portable lead- 
acid batteries have automotive use, use in motorcycles, 
railroads, all types of military vehicles on the land, sea, 
and air, and additional uses not related to vehicles. 

Additionally, there is a wealth of specific information 
concerning lead-acid batteries in literature. Volume after 
volume has been written exploring the numerous characteristics 
of this type of battery in great detail. Chapters upon 
chapters have been written about the history, electrochemistry, 
applications and performance characteristics of lead-acid 
batteries. Detailed studies have been done on performance and 
structural aspects, as well as on self-discharge effects. 

Because the author chose as the aim of his thesis, the 
establishment of an electrochemical model of the cell, based 
on Shepherd's equation (1965), it was necessary to chose a 
battery type for the experimentation which was rechargable, 
lightweight, and portable, in addition to having similar 


Characteristics of other types of cells. The lead-acid 
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battery, in addition to being rechargable, is available in a 

clear plastic case, thus affording the investigator a visual 

means to determine some battery conditions which are not 

readily available empirically in a battery with an opaque case. 
Also, it was felt that the characteristics of some of 

the other types of cells are similar to those of the lead-acid 

cell and that a model for the one may be applicable to the 


others. 


1.2 HISTORY OF LEAD-ACID BATTERIES 

The original incentive toward the development of an 
electrochemical source of current was the "galvanic cell" 
which was based on Galvani's "frog leg" experiment, circa 
1789. 

Volta, after 1792, developed several useful current 
sources, such as the Volta pile. He also formulated the 
first electrochemical theory on the subject. 

Ritter developed a cell in about 1803, Ritter's pile, 
which, when disconnected from an external source of current, 
would yield a current in the opposite direction. This was the 
basis for the secondary cell or storage battery. 

In 1860, Planté demonstrated his cell in which two lead 
Sheets separated by rubber were rolled into a spiral and 
placed in a 10% sulfuric acid solution. When charged, the 


amount of energy depended upon the amount of lead oxide formed. 
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By 1880, W.V. Siemens had invented the dynamo, which 
was then installed in central power stations. Batteries were 
used for load leveling purposes, being charged with energy 
produced at night. 

Around 1881, Faure coated the lead foil with lead oxide 
for the purpose of accelerating Planté's forming process. 

He used various lead current collector grids. 

1883 brought us the Planté plate which was actually 
designed by Tudor. In this design, he pasted lead oxide onto 
cast ribbed purelead plates. Grid plates, plates cast from 
antimonial alloys were also developed at about the same time. 

From 1890 to 1920, wood was used to separate the plates 
and then rubber. Around 1927, porous ebonite and microporous 
rubber was used as separators. 

Around 1965, batteries were designed to be relatively 
maintenance~free for use in portable items, and, with the 
advent of plastic battery cases and lightweight construction, 
automation and mass~production became feasible. 

The 1970s gave us new applications of the lead-acid 
battery, with increasing interest in Traction batteries coming 
to the forefront. Water-activated and maintenance-free 
automobile batteries were designed. Kordesch, (1977) 

These historical developments, some important, some 
merely historically interesting, represent the basis for the 


lead-acid battery as it is known today. 
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1.3 THE NEED FOR A MODEL 

A system designer needs an accurate electrochemical 
battery model in order to produce a design which makes optimum 
use of these batteries. With the use of computer-aided 
design, a functional description of various battery character- 
istics would greatly simplify the present cumbersome task of 
establishing an algorithm for interpolation of battery 
characteristics from its family of discharge curves. An 
algorithm for predicting the response to an arbitrary (non- 
constant current) discharge or charge would be, at the very 
least, extremely difficult, if not prone to significant 
inaccuracies. 

Several empirical models have been developed in the past 
to show the dependence of capacity on discharge current I 
and discharge time t {Schroder (1894), Liebenow (1897), 
Peukert (1897), Rabl (1936), Davytan (1946), and Shepherd 
(1965) }. Shepherd's (1965) model (discussed in detail in 
Chapter II) gives a relationship between cell voltage on dis- 
charge as a function of current density and time based on the 
theoretical laws of electrochemistry. Shepherd's equation 
(1965) serves as the basis for this investigation into estab- 
lishing a State Model. 

Very early in this investigation, it was found that the 
equation provided a satisfactory fit to a single discharge 


curve (complete discharge at constant current) but a poor fit 
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for the family of discharge curves. A typical family of 
discharge curves is shown in figure 1-l. Also, Shepherd 


did not include anv term which would account for secular life 


changes in the cell. 


Voltage 





0 20 4 0 6 0 B 0 100 
capacity, Ah 


Fig. 1l-l. Discharge Curves for a Typical 100-Ah (20-h rate) 


SLI Lead~Acid Battery. 


The emphasis of this thesis was to establish a model 
which would more accurately predict the discharge characteris~ 
tics of a cell at any arbitrary current, thereby providing a 
better fit for the family of discharge curves. Basing his 


investigation into a State Model on the Shepherd equation 
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(1965), the author has chosen to call his model the "Modified 
Shepherd Equation." It was also the desire of the author 

to experimentally determine the secular life changes of the 
lead-acid cell and to include a model of these changes in the 
Modified Shepherd Equation. 

The problem of needing an adequate model was made clear 
to the author while working on a Diesel submarine computer 
synthesis model at the Naval Ships Engineering Center in 
Washington, D.C. A great deal of time was invested in estab- 
lishing a model of the lead-acid cell which would adequately 
predict the capacity of the cell at arbitrary discharge cur- 
rents. This was accomplished simply as a curve fit to the 
capacities determined from the family of discharge curves. 

Due to the complexity of the problem, no attempt was made 
to use the discharge characteristics in that model. If a model 
such as the one proposed in this thesis had been available, not 
only would a significant number of man-hours have been saved, 
but also another degree of flexibility would have been allowed 
in the synthesis model. 

Also, the problem of the design of an electric vehicle is 
essentially identical in this respect. Moreover, the battery 
Manufacturer, by supplying only a few coefficients for a 
battery model would be providing the characteristics of his 
product in a very concise form which would require no inter- 


polation for the user's specific application. 
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1.4 BATTERY USAGE 

Since the earliest "galvanic cell", the lead-acid storage 
battery has been in demand for numerous purposes. Marine 
uses alone are diverse. 

Marine batteries, which are rugged and durable, have 
given long service to ships of the Navy and Merchant Marine. 
These batteries function under extremely adverse conditions 
and are subject to continual vibration. 

These batteries have been used in all types of ships from 
the largest aircraft carrier to the smallest pleasure craft to 
power radios, lights, pumps, hoists, and compressors. They 
are used in torpedoes and to fire guns. They are used in the 
kitchens of ships to provide power to stoves and refrigerators. 

Batteries power alarms and emergency lighting and even 
pump bilges and hoist anchors. Many are used in ignition 
Circuits of internal-combustion engines of small craft. They 
provide emergency power for submerged Nuclear submarines and 
the propulsion for submerged Diesel submarines. 

Marine batteries are used for gyrocompass emergency power 
and for telephones, and emergency generator start-up/ or 
field flashing. 

For many years, scientists and engineers have studied the 
possibility of converting from gasoline powered motor vehicles 
to electrically powered motor vehicles, from hybrids to all- 


electric delivery vans. Much interest was in all-electric 


o> 
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automobiles. It has been found that there are no savings in 
energy if batteries simply replace the gasoline engine. 
Electrically powered vehicles would offer a significant 
pollution advantage over gasoline and Diesel engines. 
Batteries could be used to pvrovide additional energy require- 
ments in hybrid gasoline/electric vehicles for load leveling 
purposes or for applications where the possible recovery of 
energy through regenerative breaking is significant. 
(Kordesch 1977) 

Batteries have been classified according to usage, 
according to Bode (1977): 

il SLI Batteries (Starting, Lighting and Ignition) 

These batteries are used to start an internal- 

combustion engine and to provide power to the 

electrical system when the engine is not running. 

The largest number of storage batteries are for 

use in the automobile and are used for this purpose 

more than for any other. 

te Stationary Batteries 

These generally provide emergency power when 

the primary source of power is not available. 

These include such applications as communications 

systems, electric utilities, computer systems, 

emergency lighting, and railways to provide peak 


loads, emergency power or smoothing of ac power. 





Le 


Be: Motive Power (Traction) Batteries 

Lead storage batteries provide power for the 
propulsion of devices such as electric lift trucks, 
Mining equipment, street delivery vehicles, and 
other types of material handling. 
4, Special Purpose Batteries 

Some examples of batteries in this category 
are aircraft, submarine, special military and small 
sealed batteries for consumer applications. Some of 


these uses have already been discussed. 
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me BATTERY LIFE 

The life of a battery can be defined as the period of 
time during which it has the ability to perform within certain 
specified limits. In other words, the end of life occurs 
when the battery can no longer perform within these limits. 
These limits actually depend upon the apvlication and the re- 
quirements of the system they serve. The end of life of an 
automobile battery occurs when it no longer is able to provide 
sufficient cranking power to start the engine. For a batterv 
powered vehicle, it may be when the vehicle will no longer 
travel a certain number of miles at a specified speed; in 
traction batteries, the decrease of capacity below a certain 
value: may be the criterion on which the end of life decision is 
based. Generally, however, systems are designed to meet all 
specified requirements as long as the battery is able to 
deliver 80% or more of its new capacity. 

Bode (1977) states that the end of battery life is best 
evaluated in the laboratory using fixed standards and speci- 
fied test conditions. The ideal requirement is that the life 
test correspond to the load conditions of the battery in the 
field, however, as in the case of automobile batteries, ideal 
testing in the laboratory is impossible because it is im- 
possible to duplicate such diverse conditions as automobile 
type, driving methods, maintenance and environment. 


Pigure 1-2 is a plot of the change in capacity during the 





24 


lifetime of a lead-acid battery. Very little information is 
available which provides details on the change of characteris- 
tics of a cell during its lifetime. It is therefore proposed 


that this data be obtained by experimental means. 






2D Op ao ee a=) 


% Of Nominal Capacity 





oS ee eee ee ee eee ee ee ee 


vei en pe—neme 


End of 
eee es Life Point 


Figure 1-2. Typical Cell Capacity During the Lifetime of 


a Lead-Acid Battery. KORDESCH (1977) 


The discharge capacity of a battery has been found to be 
somewhat dependent upon the rate of previous discharges. Bode 
(1977) calls this the "memory effect." Vinal (1955) states 
that “the capacity is lower if the discharge has been preceded 
by a discharge at a higher rate, and it is higher if preceded 
by a discharge at a lower rate." Therefore, any experimental © 


procedures should take the "memory effect" into consideration. 
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1.6 REASONS FOR DECREASE IN LIFE EXPECTANCY OR FAILURE 

According to Vinal (1955), storage batteries large and 
small and in widely diverse services are subject to conditions 
that may lead to excellent service or through accident or 
abuse to something less than satisfactory. Under his heading 
"Sources of Trouble," he writes about faulty conditions 
which if detected could be corrected. Detection and correction 
of faulty conditions might prolong the satisfactory battery 
life into years. 

He also states that as a result of studies done in 1947 
on thousands of worn out batteries which had been used in 
normal service, the reasons they wear out became well known. 
In comprehensive data relating primarily to automobile 
batteries, grid corrosion of the positive plates accounted for 
42% of the total failures. Naturally, with the use of new 
alloys, batteries developed since that time should have less 
positive plate corrosion. 

Bode (1977) indicates that destruction of the positive 
grids occur when lead sulfate forms a heavy layer on the grid 
to protect it from being dissolved anodically. One rarely 
encounters corrosion of the negative grid. Hydrogen corrosion 
occurs in the presence of organic substances and at higher 
Operating temperatures. The result of corrosion is loss of 


conductivity and grid breakage. 
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Vinal (1955) indicated that 7 to 10% had less than total 
life span due to shedding of positive active material, 
cracked partitions and leaking cases accounted for nearly one- 
third, and buckled plates and short-circuited separators pro- 
duced failure in 10 to 12% of the batteries. About 10% of the 
failures were caused by negative plate failure. Sulfation 
which causes buckled plates and negative failures, was a cause 
of battery failure in less than 20% of the batteries. 

Overcharging loosens active material, corrodes positive 
grids, causes excessive gassing and produces sludging of a 
fine brown sediment. It may cause raised temperatures contri- 
buting to battery destruction with damage to plates and separ- 
ators, buckling of plates and needless loss of water. 

Bode (1977) states that in sulfation, - "hard" lead 
sulfate forms if chargeable lead sulfate undergoes long or high 
temperature aging process, (recrystallization). Sulfation 
produces loss of capacity and is caused by long standing in a 
discharge condition, too high acid concentration, increased 
self-discharge, continuous operation between 45° and Be, 
and prolonged undercharging. 

He also indicates that sludging of the positive mass 
occurs most at the end of a charge and at the beginning of a 
discharge. Heavy gassing facilitates sludging. There is a 
Qreater amount of sludging in discharges with large current 


density. The type of discharge influences sludging. If acid 
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concentration is drastically reduced, temperature increases or 
current density is diminished, there is less sludging and 
the battery life is increased. 

Vinal (1955) concludes that occaSional overcharging is 
beneficial, but habitual overcharging decreases useful battery 
life. 

He indicates that undercharging gradually runs down cells 
producing diminishing values of specific gravity. Plates be- 
come lighter in color and a sludge, a fine white powder which 
is predominantly lead sulfate, is produced. One or more cells 
could become exhausted. He concludes that the remedy is to 
charge until all cells are in normal condition. This is refer- 
red to as an equalizing charge. Undercharging is the most 
common cause of buckling of plates because lead sulfate occu- 
pies more space than the original material producing excessive 
Strain on the plates. 

Failure of the charging system and of the battery ulti- 
mately is blamed on the battery. When an automobile reaches 
speeds greater than or equal to 12 to 15 mph the battery begins 
to charge. If the generator doesn't come on at the proper 
speed the battery won't receive its required amount of charge. 
Also, the system must regulate itself properly so that the 
battery will not receive massive amounts of current when the 


automobile reaches highway speeds. Conditions in the charging 
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system which drain the battery and cause it to remain dis- 
charged are grounded circuitry and partial grounds. 
(Vinal, 1955) 

Other possible causes of shortened battery life include 
corroded terminals, short circuits, separator failures, 
worn-out plates, electrolyte below the tops of the plates, 
freezing, when the battery is in a discharge condition, 
impurities in the solution, spalling of the active material, 
growth of positive plates, shrinkage of negative plates, 
leading, and battery explosions. The interested reader is 
directed to the classic text of Vinal (1955) and the excellent 


texts of Bode (1977) and Kordesch (1977). 
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CHAPTER 2 
THEORETICAL ANALYSIS 
2.1 ELECTROCHEMISTRY OF THE LEAD-ACID BATTERY 
The theory of operation of the lead-acid battery is 
based on the “double sulfate theory" proposed by Gladstone 
and Tribe, (1881-1883). 


discharge 


PbO, + Pb + 2H2SO, 2PbSO, + 2H20 (221) 


charge 


The principle was so named because the discharge product 
of both the lead anode electrochemical oxidation reaction and 
the lead dioxide cathode electrochemical reduction reaction is 
lead sulfate. 

This early equation reveals the stoichiometric relation- 
ship of reactants and products, as well as the role played by 
the sulfuric acid. 

A more complete representation of the actual electro- 


chemical reactions is: 


ANODE: Pb + Hso, C28¢28r9e ppso, + ot + @ (22) 
aa, 
charge 
+ -_ —_ 
METHODE: PbO, + 3H + HSO, + 2e —__»PbSO, + 2H-20 (223) 


BUERALL: Pb + PhO, + 2H’ + 2HSO, ———»2PbSO, + 2H.O (2.4) 
es 
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tee clecenolytic dissociation of sulfuric acid into 
H’ and HSO, is closer to the true picture in the cell. The 
dissociation into 2H” and SO, occurs in solutions much more 
dilute than those found in a lead-acid cell. 

During overcharge, the electrolysis of water occurs, 


yielding hydrogen and exygen. 


Based on the free-energy change of 


AG = -13.95 kcal/mole 


for the anodic reaction 


(E = 0.303 V versus H,/H’) 


and 


AG 


-74.99 kcal/mole 
for the cathodic reaction 


(E = -1,.627 V versus Hey ae 
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Phe overall reaction values are: 


AG = ~88.94 kcal/mole 

(E, = 1.93 V) 

specific energy = 161 Wh/kg 
specific capacity = 83.5 Ah/kg 


Bode (1977). 


In practice, significantly lower specific energies are 
exhibited, but the theoretical values provide an upper limit to 
possible future developments. For a complete discussion of 
electro-chemical theory, the reader is directed to the 
excellent text by Bode (1977), and the classic text in the area 


By Vinal (1955). 
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2.2 THE SHEPHERD EQUATION 

Shepherd (1965) has proposed a static model for cell 
discharge characteristics which uses simple analytical 
expressions of chemicophysical origin and describes the cell 


potential E as follows: 


E(q,i) = Eo - Vv, ¢a) SV gee) oS Vq (i,q) (25) 


where Eg= open circuit cell potential 





Vv, (a) = A exp (- ) = activation volarization 

Voli) = Roi = ohmic polarization 

Vq4(i,q) = K ( 5 ) i = diffusion (concentration) 
polarization 

i = the discharge current 

q = it = the accumulated charge 


the cell capacity 


K 


the mean polarization resistance 


Ro = the internal cell resistance 


Shepherd used "it" in place of "gq" in his equation 
and Ro corresponds to his N. The activation polarization 
term is the single term which does not have any particular 
theoretical basis. Its inclusion accounts for the rapid 
drop in potential at the beginning of discharge. Figure 2-1 


shows an exagerated view of this phenomenon. The time 
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duration of this initial drop was found to be in order of se- 
conds. Consequently, the description of cell discharge should 
not be significantly affected by simply ignoring the term; this 
was the approach used by Shepherd, also. 

Another phenomenon shown in fig. 2-l, not treated by Shep- 
herd, is known as the "voltage dip" or “coup de fouet" where 
initially the voltage drops below the "steady state" value and 


then rises to this level within a few minutes. 


Voltage 


Ee de fouet! 





Time 


Figure 2-1. Cell Voltage versus Time for a Constant 
Current Discharge 
The ohmic polarization term accounts for the ohmic resis- 
tance through both the metal and the solution. This has been 
modeled as a constant, however, it is felt that this assump- 


tion may be errcneous. 


\ 
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The specific electrical resistance of lead oxide is on 
ene order of 107° Q/m and metallic lead, on the order of 
10~’ Q/m. The electrolyte has a range of specific resistance 
of about 125 2/m when nearly fully charged to about 170 2/m 
when discharged. The product of the two reactions, lead 
sulfate, is highly resistive, specific resistance on the 
order of 10!° Q/m. Based on these values, it would be expect- 
ed that the ohmic resistance increase as the reaction proceeds 
and becomes fairly significant as it approaches completion. 
The diffusion (concentration) polarization term appears 
to have been made up from a combination of what Bode (1977) 
calls the charge transfer overvoltage, Nee and the concen- 
tration overvoltage, Nee Following Bode's development, 


the terminal voltage of the cell is 


E = &£E 7 Gi 2 at) aa Pa 


rev Cun t Nolo 7 IR (2.6) 


De ay is the open circuit voltage and IR includes the 
resistances of the electrolyte and separator. The subscripts 
A and C refer to the anode and cathode overvoltages, respec- 
tively. All of the overvoltages decrease the cell voltage on 
discharge, and overvoltages of the electrodes are related to 


the current densities, 1. Assuming that the electrodes have 


equal surface area Ay = Aa = A, then for both electrodes the 





a2 


following expressions apply: 





be = Za og iy (6) b log i (ZeF) 
ZF 
and 
_ 23" RE log (1 _ 1 ) . (2.8) 
Jouay 25 1 


For small current densities, these become 





Ne — = 5 See (2 9) 
ZFig (al -Q) 
and 
_ RT i 
lek= op “Hip, (2.10) 


For this case, one of the terms in (2.6) may be neglect- 
ed if they differ by an order of magnitude. The terms in the 
equations above are: 

1) = the apparent exchange density 


the number of electrons transferred in the 


Z 


determining step 


ry 
Il 


the Faraday (measure of capacity) 


oe) 
Il 


the gas constant 
T = the temperature 
1, = the diffusion limiting current 


Z 
and 


© = the degree of coverage, (only the uncovered portion 


(1-0) of the electrode can react) 
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Replacing the constant terms in Ny by a single constant k 
yields 


— K (221) 


(1. =) 

If the uncovered portion of the electrode is related to 
the cell capacity Q and accumulated charge q where q < Q 
with gq > Q at complete discharge, equation (2.11) becomes 


Meee Ae. (2.12) 


a 
: (1-0) (Q-q) 


This term is similar to Shepherd's diffusion polariza- 
tion but with the striking difference that it is not a 
function of current. If one were to fit Shepherd's equation 
to a single discharge curve, this difference would not be 
apparent because the product of K and i in equation (2.5) 
would be equal to K in equation (2.12). 

Prior to the discovery of the above relationship, the 
author, while attempting to make a fit for a family of dis- 
charge curves to Shepherd's equation, produced fits of each 
curve individually. In the process of analyzing these results 
a rather surprising discovery was made in that the product 
of K and i was nearly constant. 

Another item in Shepherd's equation which should be 


addressed, is the cell capacity, Q. The value for Q should be 


close to the capacity of the cell at the particular discharge 
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current. In fact, it needs to be slightly larger than the 
actual capacity to avoid predicting infinite overvoltages, 
i.e. Q > q for all values of g at the particular discharge 
murrent. 

The most widely used model for all capacity is Peukert's 


(1897) equation, 


t=c. (22S) 


In this equation n and C are constants which can be determined 
by using the discharge rate and time for complete discharge 
at two different currents. Since the cell capacity in terms 


of Q is the desired quantity, the substitution of t = Q/I gives 


pM, =C or O= cy finn) (2) 


The values of n and C can be calculated as follows: 


(Sol) 


OTea oo, 

Q,=¢ Eye 

(n-1) log tee ee arog On 
(n-1) log Og C= og O- 


TOqsOsm = 1090, 
n= ——————————- + ]l 
Weg =, — log tf, 
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The value for n is in the range of from 1 to 2. For 
diminishingly small current densities, n approaches 1 and 
for high current densities, n approaches 2. Typically, the 
value of n falls within about 1.3 to 1.4 This equation is 
not based on any particular physical laws, but serves as a 
good empirical fit to discharge capacities. 

Baikie, Gillibrand, and Peters (1972) expanded on this 
relationship to include temperature effects in the range of 
-18°C to +16°C; for current densities between 21 and 


165mA/cm* it has the form 


Where Q is the effective plate capacity, and &@ is the temp- 


erature in . 


C. The values they established were 
n= 1.4 for both plates 
Kee Oes2 and 0224 min, 


and a = 0.021 and 0.015 min/°c 


for the negative and positive plates, respectively. 

Even if these equations provide a perfect correlation to 
the capacity of a cell, they do so only at one instant and 
do not predict changes which occur throughout the life of the 
battery. If a state model is to be established for the electro- 


Biemical cell, then these effects must also be included. 
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eS METHOD OF DETERMINING COEFFICIENTS FOR THE MODEL 
An iterative approach using a digital computer was used to 
determine the coefficients for the model. This curve fitting 


method basically used the following steps: 


18 Choose initial values 
Die Vary each value independently 
Trial Value = Value * nee a 
where n = -1,0,1l and 0O Factor 1 
oe Choose the set of trial values which gives the 


smallest difference between the data and model. 


Steps 2 and 3 were repeated until no further improvement 
was required by the user. When the "best" set of values was 
found, then these trial value were used in the next iteration 
as the values. When a particular factor failed to make any 
improvement, the difference between the factor and 1 was reduc- 
ed by 1/2 and a new iteration was performed. 

For runs of the Shepherd equation, initial values were 


chosen such that: 


E, = Open circuit voltage 


Q = cell capacity plus a small amount to prevent 
division by zero 


Pele Rt = EQ — discharge voltage at t = 0. 
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This method finds a minimum value for the correlation, 
but the only assurance is that the set of values corresponds 


£O a relative minimum. 
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CHAPTER 3 
EXPERIMENTAL PROCEDURE 
3.1 BATTERY CHARACTERISTICS 

In order to establish a model which would more accurately 
predict the discharge characteristics of a cell at any arbi- 
trary current and throughout its life, testing was verformed 
using a lead-acid battery. 

The type of battery used was a low capacity motorcycle 
type battery in a clear case. The low capacity was used so 
that the equipment could be chosen requiring smaller current 
carrying capacities. The advantage of the clear case over an 
Opaque case was that the author could more readily observe 
conditions of.the battery, including the amount of gassing 
during charging, the level of the electrolyte, the changes, if 
any, in the condition of the plates and the accumulation of 
any sludge, or fine-grained sediment, in the bottom of the bat- 
tery case without destroying the battery in order to make these 
observations. 

Table 3-1 lists the battery specifications and overall 
dimensions. Figure 3-l shows the internal configuration of a 
Cell. Each cell was configured with anode (-) plates on the 
ends: As can be seen, in the figure, there are four anode (-) 
plates and three cathode (+) plates with 6 separators in the 


Spaces between the plates. 
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Internal Configuration of a Cell. 


Figure 3-1. 


Wisco Division 


Manufactured for 


WI 53403 


ESB Incorporated 
1222 18th Street 


Racine, 
JA-11 


3 Cell 


6V 


: 
; 
i 


E6N6—-3B 


JIS Designation: 


-h rate 


6Ah @ 10 


Clear Plastic Bottom 


with Black Plastic Ton 


External Case Dimensions: 


Battery Specifications and Dimensions. 


Table 3-1. 
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Connector 


Note: All plates 

appear to have the same 
Anode Plate #1 dimensions. 

Cell #1 Battery #1 


Area x 29cm 


62 
TAM 


Note: This tab suv- 
ports the plate ona 
rib in the case run- 
ning the length of the 





Figure 3-2. Side View of Anode Plate #1 as Seen from the 


End of the Battery. 


Table 3-2 lists the measurements on cell #41 to de- 
termine the thicknesses of the vlates and separators. These 
measurements were made by placing a rule against the side of 
the case. Some degree of uncertainty in the measurement of the 
plate thicknesses can be attributed to the method of measurement 
required. 

Figure 3-2 shows the side view of anode plate #1 which is 
the view seen when looking at the end of the battery. All 
plates have essentially the same dimensions. The area of the 


plates is approximately 29cm’. 
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COMPONENT DISTANCE THICKNESS 
(rm ) (mm) 

Anode Plate #l TS 1 eas. AVERAGE THICKNESS 

Separator #1 4.0 2.5 Anode 157 Sm 

Cathode Plate #1 Se2 Dae S. Dev. Oe222 

Separator #2 oO 2.8 Var. 0.0369 

Anode Plate #2 10.3 divers 

Separator #3 35 0 247. Cathode 2.07mm 

Cathode Plate #2 15.0 2.0 S. Dev. OulaS 

Peparator # 4 NT od 7G) " Var. 0.0089 

Anode Plate #3 195 12 

Separator #5 220 225 Separators 2.62mm 
_ {Cathode Plate #3 24.0 Pere) S. Dev. eke 3 

Separator #6 26.5 23 Var. Ones 

Anode Plate #4 232 te 7 


Table 3-2. Measurements of the Distances of the Components 
from Edge of Anode Plate #1 and the Thicknesses 


of the Anode and Cathode Plates and Separators. 
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Figure 3-3. Top View of Test Battery 


Table 3-3 shows the volume of the materials in a cell, 
including the deduced volume of electrolyte of 80cc. 

Access to the connectors between cells #1 and #2, and 
cells #3 and #4 was made by drilling into the tov of the case. 
The location of the holes is shown in figure 3-3. 

The following observations were made of the battery in its 
unused condition. The anode plate viewed from the side had a 
light grey appearance; viewed from the end it was darker grey. 
Bubbles were seen on the plate and periodically several would 
rise to the surface. Cathode plates could only be observed 
from their ends due to the configuration of the cells. They had 


a dark brown rusty coating on their surfaces. Bubbles were not 
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INTERNAL VOLUME (per cell): 


L = 29mm 
W 51.5mm 


Height to upper electrolyte level (marked on case) = 80mm 


Volume = 120cc 


VOLUME OF ELECTROLYTE: 
Vol. of anode (per cell) = 29cm « .1575cm *« 4 = 18.3cc 
Vol. of cathode (per cell) = 29cm * .207cm * 3 = 18.0cc 


Internal volume - plate volume = 120cc - 36cc = 84cc 


Assume solid volume of separators = 4cc 


VOLUME OF ELECTROLYTE = 80cc/cell 





Table 3-3. Volume of Materials within a Cell. 
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generally observed on these plates. 

These observations were made several weeks after the 
batteryhad first been activated, but before it had been 
used. When the battery had first been activated, the anode 
and cathode plates had more nearly the same appearance, which 
was light grey. 

The separators appeared to be made of a "fiber wool” 
placed on a paper (or felt) backing. The backing was always 
placed against the anode plates, while the "wool" was against 
the cathode plates. Small bubbles were entrapped in the fiber. 

The cavity below the plates had a very slight amount of 


rusty-colored, fine-grained sediment. 


3.2 EXPERIMENTAL SET-UP 

Figure 3-4 is a circuit diagram of the experimental set-up 
using a double pole-double throw switch. The circuit was 
configured in this manner so that no disconnecting and 
reconnecting would be required between charges and discharges 
and a single ammeter could be used in either mode. 

Two different battery charging set-ups were used. One was 
simply the use of a commercially available 6 Volt Tapering-Type 
battery charger. The battery charger was manufactured by 
Mallory, Indianapolis, Indiana (Engineering Projects Lab 
Number 1167), Model 6SAC1O. 


Two chart recorders were used. One was a 4 channel chart 
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Figure 3-4. Circuit Diagram of the Experimental Set-Uo 


Using a Double Pole-Double Throw Switch. 


recorder manufactured by Gould Brush 2400 Model 2007-4404-00 
Serial 00473. All four channels were used to record the 
voltage on each cell and the current simultaneously. This 
provided the experimenter with a complete picture of the 
behavior of the battery. 

The single channel chart recorder, a Hewlett-Packard 
Mosley 680 Strip Chart Recorder, Model H01-680, Serial # 
7112-01825 (Engineering Projects Lab #1175) was used to record 


the voltage of the cell with the lowest capacity, i.e., the 
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cell which controlled the end of discharge point. This was 
cell #3 in battery #2. The single channel chart recorder was 
used twenty-four hours a day during life cycle testing, 
primarily because its slowest speed of 1 inch per hour gave a 
concise record of charge/discharge cycling. 

The four channel chart recorder was used when a more 
complete picture of the charge/discharge cycles was desired. 
The slowest speed on this chart recorder was 5mm ver 100 sec 
(approximately 7 inches per hour). Continuous running would 
have produced an unmanageable amount of paper. 

The set-up for life cycle testing included a sensing 
circuit to determine the end of discharge and end of charge and 
a relay to perform the switching of the battery circuit, as 
shown in figure 3-4. Figure 3-5 shows the schematic of the 
sensing and relay-driving circuit. The sensing circuit is a 
schmitt Trigger described by Millman and Taub (1965). 

It is a bistable circuit, and for a loop gain greater 
than 1, its input-versus-output displays hysteresis, as 
shown in figure 3-6. 

The relay driving circuit simply consists of a single 
transistor connected in series with the relay, with the out- 
put of the Schmitt Trigger to the base of transistor Q3. The 
diode connected across the relay was simply for the purpose of 


providing positive opening of the relay when the driving 
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Figure 3-5. Schematic Diagram of the Schmitt Trigger 


Sensing and Relay-Driving Circuit. 


transistor turns off. The other relay shown in the circuit was 
included as a fail-safe device. It is energized by 115V ac 
line power. Its terminals are connected in series with the 
ammeter so that if there is a loss of power, the battery cir- 
cuit is opened. The terminals of the main relay were wired to 
the terminals of the double pole-double throw switch through a 
plug. When it was desired to perform controlled testing of 
the battery, the sensing and relay circuit was unplugged. 

One recommended modification to the circuit is to insert a 
B position switch to the base of the driving transistor Q3. 


One position would connect the Schmitt Trigger to the tran- 
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Figure 3-6. Hysteresis of Schmitt Trigger Circuit. 


Sistor, as shown in figure 3-5. Another would bias Q3 
off, placing the hattery on discharge, while the third pos- 
ition would turn Q3 on placing the battery on SEES. This 
would facilitate making the controlled charges and discharaes. 
Also, a simple open or closed switch in the fail-safe relay 
line might be added. With the set-up as shown, it was 
necessary to unplug the relay's energizing power from the wall 
receptacle. 

The Schmitt Trigger was very sensitive to fluctuations of 
Short time duration. This might be improved by inserting small 


Capacitors from the base of Ql and Q2 to their emitters. 
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3.3 TESTING PROCEDURE 

There were basically two objectives of the testing. 

One was to experimentally obtain a family of discharge curves 
for the verification or modification of the Shevherd (1965) 
equation. The other was life-cycle testing in order to 
determine the change in battery characteristics with age. 

To obtain a family of discharge curves, constant current 
discharges were performed at 0.6, 1.5, 3.6, and 5.4A. These 
values were chosen because it was expected that they would 
provide a range of values similar to those shown in a plot 
of typical curves, such as those in figure 1l-l. Once two runs 
have been made at different currents, Peukert's (1897) 
equation can be used to determine the eee recuired for the 
discharge times desired. 

In an effort to place the battery in identical conditions 
prior to a discharge in which the data would be used for 
comparisons, the following procedure was followed. The 
battery was fully charged, followed by a pre-conditioning 
discharge at 3.6A. The reason for the pre-conditioning dis- 
Charge was to remove variations in the "memory effect" (dis- 
cussed in Section 1.5) caused by previous discharges. The 
pre-conditioning discharge was followed by a constant current 


Charge of 3.0A. The constant current was maintained by the use 


© 





oye 
Voltage 
3.0 


Time 


Figure 3-7. Typical Charging Curve at Constant Current 
Of 20k. 

of a 12 Volt battery connected in series with a variable 

resistor and was used as a charger. 

These values were chosen because they were within the 
range of currents used during life cycle testing and for con- 
venience. By using this scheme, it was felt that repeatable 
results could be obtained. The end of charge was determined to 
occur when the battery was gassing heavily and the cell voltage 
had remained constant for about 10 to 15 minutes, as shown in 
figure 3-7. 


Prior to discharging the battery, it was left on open 
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circuit for a period of at least one hour but not generally 
more than twelve hours. This would allow the electrolyte to 
become fairly uniformly distributed throughout the cell and 
also for the temperature to stabilize. 

During discharge, the variable resistor which served as 
the load was adjusted periodically to maintain constant 
current. The end of discharge occurred when the voltage of a 
cell dropped off abruptly. Discharge was terminated when this 
voltage dropped to approximately 1.0V. 

Life cycle testing was performed using the Schmitt 
Trigger circuit and relay with the Mallory battery charger 
and a constant load resistance of about 1.62. The charging 
Burrent was in the vicinity of 3.0A, ranging from 3.5 at the 
start of charge to about 2A at the end. The average discharge 
current was in the vicinity of 3.5A. The Schmitt Trigger 
circuit was set to end discharge and commence charging when the 
battery voltage dropped to about 5V, and to switch on the high 
end at about 7.5V. The setting at the low end could vary by as 
much as several tenths of a volt without any particular change 
an the time of switching, due to the very rapid fall off of 
voltage at this point. The setting at the high end needed to 
be much more precise due to the leveling off of voltage at the 
fully charged point. This proved to be most difficult and 
introduced variations in the state of charge at the commence- 


Ment of discharge under life cycle testing. 
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CHAPTER 4 
RESULTS AND DISCUSSION 

4.1 MODIFICATIONS TO THE SHEPHERD EQUATION 

The numeric data corresponding to figures 4-1 through 4-36 
in this section is found in appendix A. Figure 4-1 shows the 
family of discharge curves obtained following the methods out- 
lined in Section 3.3. Figure 4-2 shows the same curves after 
manual smoothing. One of the primary effects of this smoothing 
was the removal of the initial rapid drop of voltage, which 
Shepherd accounted for in his activation polarization term, and 
the "coup de fouet." Since the initial drop was completed 
within seconds of starting the discharge, and the "voltage dip" 
had disappeared within about two minutes, ignoring these pheno- 
mena should have little effect on the applicability of a model. 

The first step toward the establishment of a model was an 


attempt to fit the discharge data to Shepherd's equation: 





Q ; 
E = —s —_ 
E. K O- q a Rol 


The values found by a least squares fit were: 


Bee 2295 
S 

K = 0.08086 

O= 6.844 


Ro = 0.00092 


with a correlation © (E - eo) for all curves = 3.5008. 
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Figure 4-1. Family of Experimentally obtained Discharge 
Curves. 
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Figure 4-2. Smoothed Family of Discharge Curves. 
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Figure 4-3 shows the resulting family of discharge 
curves. Comparisons of each individual discharge curve and 
its corresponding modeled discharge curve are shown in figures 
4-4 through 4-7. 

These results illustrate the need for some modification. 
The rapid drop-off of voltage at the end of discharge was only 
predicted for the .6 amp discharge. This was the peSEeA 
result, considering the denominator of the concentration 


polarization term. 
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Figure 4-3. Family of Curves for Curve Fit to Shepherd's 
Equation. 
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Figure 4-4. 0.6 Amp Discharge Curve and the Modeled Curve for 
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the Fit of the Family of Gurves to Shepherd's Equation. 
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Figure 4-5. 1.5 Amp Discharge Curve and the Modeled Curve for 


the Fit of the Family of Curves to Shepherd's Equation. 
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3.6 Amp Discharge Curve and the Modeled Curve for 
of the Family of Curves to Shepherd's Equation. 
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5.4 Amp Discharge Curve and the Modeled Curve for 
of the Family of Curves to Shepherd's Equation. 
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Next, a curve fit for each individual curve was per- 


formed using the Shepherd equation. 


Table 4-1 lists these 


coefficients with the graphical results shown in figures 4-8 


through 4-1l. 





Current, i .6 eo 3.6 Sve 
E. 18 alles Pape Us Ze & 
K Oe Des -02358 0.0031 0.00142 
0 yo02 2202 ec SMe Paco Al 
Ro PEI SLS, .0481 0.0734 0.05924 
Correlation .0718 e239 Ozer 5 35 0.2638 
Bele 4-t. Curve Fit of Shepherd's Equation to Individual 
Discharge Curves. 
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Figure 4-8. 0.6 Amp Discharge Curve and the Modeled Curve for 


a Fit to the Shepherd Equation of this curve alone. 
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Figure 4-9. 1.5 Amp Discharge Curve and the Modeled Curve for 


a Fit to the Shepherd Equation of this Curve Alone. 
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Figure 4-10. 3.6 Amp Discharge Curve and the Modeled Curve for 
a Fit to the Shepherd Equation of this Curve Alone. 
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Figure 4-11. 5.4 Amp Discharge Curve and the Modeled Curve for 


a Fit to the Shepherd Equation of this Curve Alone. 
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If the slope of the model were adjusted then cae cor— 
relation between discharge and model would be very good. This 
suggests modifying the resistance polarization term such that 
Vo = Vol(q,i) = (agqtb)i. Ry = agqtb takes into account that the 
ohmic resistance of the cell changes as the discharge reaction 
proceeds toward completion. This model simply predicts a 
linear change in resistance. It may be found Bae is also 
a function of current since the PbSO, build-up occurs in a 
different manner at different discharge rates. 

The value for Q follows the Peukert equation very well 
when C and N have the values 5.803 and 1.2227 respectively. 
Figure 4-12 shows a log-log plot of Q versus I for Peukert's 


equation. 
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Figure 4-12. Log-Log Plot of Peukert's Equation 
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Modifying the Shepherd Equation by the inclusion of 
Peukert's capacity model yielded a fit which predicted very 
well the drop in voltage at the end of discharge. The 
results obtained did not produce a very good fit to the left 
of this point, however. 


These plots are shown in figures 4-13 through 4-17. 
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Figure 4-13. Family of Curves for Modification of Capacity 


Term. 
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Figure 4-14. 0.6 Amp Fit with Capacity Term Modified. 
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Figure 4-15. 1.5 Amp Fit with Canvacity Term Modified. 
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Figure 4-16. 3.6 Amp Fit with Capacity Term Modified. 
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Bigure 4-17. 5.4 Amp Fit with Capacity Term Modified. 
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As mentioned in Section 2.2, the product of K and i, 
found from the curve fits figures 4-8 through 4-l1, tended to 


be constant, as shown in table 4-2. 


1 O76 Us, ces: 5.4 
K O.053s7 O0236 02.003 50 0.00142 
Kei O2Ges2Z O03 57 OF OE: 0200757 


Table 4-2. Comparison of K and the Product Kei. 


This suggests that in equation 2.5, Va x Vqalari), but rather 





= = QO : ; 
Va =V 5 (q) = K O- ¢" By including both this new Va and 


capacity terms, the plots of figures 4-18 through 4-22 were 
produced. These show a definite improvement over previous fits 


for the family of discharge curves. The values determined were 


E. = 2.002 
S 

K = 0.0039 

Ro = 0.03 


with a correlation of 1.8309. 
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Figure 4-18a.. Experimentally Determined Family of Curves 


Modification of vA = Vg (a) and Cavacity are shown below: 
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Figure 4-18b. Family of Curves for Modification of 


V 


= Vg (a) and Capacity. 
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Figure 4-19. 0.6 Amp Fit with Modification ene ite Va (a) 


and Capacity. 
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Figure 4-20. 1.5 Amp Fit with Modification of V, = Vg (a) 
and Capacity. 
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Figure 4-21. 3.6 Amp Fit with Modification of V, = V,(q) 


and Capacity. 
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Figure 4-22. 5.4 Amp Fit with Modification of V4, = V,(q) 
and Capacity. 
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Each plot of figures 4-18 through 4-22 suggests that if 
the ohmic resistance were modeled as a linear function of 
accumulated charge, a much better fit could be obtained. 

Figures 4-23 through 4-26 show the results of individual 


fits for each discharge curve by modeling 


Ro = Ro (q) = RW + Ry, 


Table 4-3 lists these results. 





1 OLS eS 3706 Dad 

Ee Pi Pails sii0) 22230 Zoe o Payal 3 0 

K 0.00744 0.03453 0.00876 O500615 
R, OU ets 0.00073 0.01881 OO a7 7A 
Ry OOS 3867 0.05424 003253 0.02016 
Correlation 0705279 Oe 2685 0.03129 ONTO ei sh7/ S 
Figure 4-23 4-24 4=25 4-26 
Table 4-3. Curve Fit for Individual Discharge Curves with 


Modifications to QO, V and Ro. 
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Figure 4-23. Fit of 0.6 Amp Alone for Modification of 


Ro and Capacity. 
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Figure 4-24. Fit of 1.5 Amp Alone for Modification of 
Ro and Capacity. 
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Figure 4-25. Fit of 3.6 Amp Alone for Modification of 


Ro and Capacity. 
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Figure 4-26. Fit of 5.4 Amp Alone for Modification of 
Ro and Capacity. 
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Twe more curve fits to the family of discharge curves were 
performed using the modified ohmic resistance term and 
Peukert's capacity model. 
The first was performed using Shepherd's Va term with the 
results shown in figures 4-27 through 4-31. The second fit 
was made using the modified Va term. These results are shown 


in figures 4-32 through 4-36. 
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Figure 4-27a. Experimentally Determined Family of Curves 


Modification of Capacity and Ro = Ro(q) are shown 


below: 
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Figure 4-27b. Family CEmCcurvecewiEn Modification of Capacity 
and Ry = Ry(q). 
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Begure 4-28 0.6 Amp Fit for Modificaticn of R» and Capacity. 
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Figure 4-29. 1.5 Amp Fit for Modification of Ry and Capacity. 
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Figure 4-30. 3.6 Amp Fit for Modification of Ro and Capacity. 
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Figure 4-31. 5.4 Amp Fit for Modification of Ry and Capacity. 
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Figure 4-32a. Experimentally Determined Family of Curves 


Modifications to Var Ro, and Capacitv are shown below: 
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Figure 4-32b. Family of Curves with Modifications to 
Var Ro, and Capacity. 
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Figure 4-33. 0.6 Amp Curve with Modifications~-to Var Ro, 


ane Capac cy . 
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Figure 4-34. 1.5 Amp Curve with Modifications to Var Ro, 


and Capacity. 
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Figure 4-35. 3.6 Amp Curve with Modifications to V,, Ro, 


and Capacity. 
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Figure 4-36. 5.4 Amp Curve with Modifications to Var Ros 
and Capacity. 
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mee LLIFE CYCLE TESTING 

As was previously mentioned in Section 1.5, according to 
Bode (1977), the ideal requirement for life testing is that the 
life test correspond to the load conditions of the battery in 
the field. However, since ideal testing is impossible due to 
conditions impossible to duplicate in the laboratory, one 
must choose another type of test. Kordesch (1977) states that 
the Cycling Life Test is a test principally to evaluate battery 
internal components. It does not attempt to simulate car 
service conditions which are extremely variable. Naturally, 
this is the type of test chosen for testing the lead-acid 
battery which was used for experimentation for this thesis. 

The life-cycle testing procedure used to determine the 
Change in battery characteristics with age is outlined in 
section 3.3. Some of the pertinent data concerning life- 


cycle testing follows: 


Number of Cycles = 135 

Time on Charge = 193 hours 

Time on Discharge = 114 hours 

Total number of Ampere Hours Discharged = 375 


Amount of Water Lost Due to Hydrolysis = 12ml/cell 
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To determine changes in cell capacity during the life 
of the battery several controlled discharges at 3.6 amps were 
made at different times during the testing. Table 4-4 
shows these results. The final capacity was 94% of the pre- 
vious capacities. Further testing is required to determine 


the rate at which the capacity decreases. 


Number of Prior Cycles 25 Galt i235 
Discharge Time TZ 72 oy ares. 


Amphere-hour Capacity £.392. 4232 °:4.05 


Table 4-4. Capacity Change of Cell #3 During the Life 


of the Battery. 


A great deal of data was generated during the testing but 
it was considered far to voluminous to be included in its raw 
form in this report. This data, which is primarily in the form 
of strip chart recordings, was placed on file with Professor 


H.M. Paynter at the Massachusetts Institute of Technology. 


4.3 UNCERTAINTIES IN EXPERIMENTAL RESULTS 
Uncertainties involved in the experimental results were 
due partly to imperfect knowledge of the state of charge and 


the inability, at least with the experimental set-up used, to 
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control certain parameters. The most notable of these was 
temperature. Variations in temperature have been recognized 
to produce significant changes in battery characteristics. 
Great variations in temperature were not encountered since the 
experiment was performed at room temperature, but there were 
some variations in this temperature by nearly as much as KO) we 
The effect of internal heating at higher discharge rates may 
have affected results. 

The difficulty of knowing the precise state of charge also 
added a degree of uncertainty. Included in this was the 
"memory effect," since the amount of preconditioning required 
to eliminate it was not known. 

Also, the knowledge of the mechanism of lead sulfate 
build-up at the various discharge rates was not available. 

Due to the differences of the results during curve fitting 
for the 1.5 amp discharge with respect to the other curves, 
there was some suspicion that the conditions for that particu- 


lar discharge may have been different than for the others. 


4.4 ADDITIONAL OBSERVATIONS 

The method of charging the battery and the point at which 
cnarging was terminated was deemed to be the most important 
factor in placing the battery in the same state-of-charge 
before the controlled discharges. Figure 4-37 shows a typical 


constant current charging curve for this experiment. 
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Figure 4-37. Typical Constant Current Charge. 


The amount of gassing in the battery was determined by 
visually observing the amount of gas evolution. The clear 
battery case was vital in making this determination. By 
observing the top of the plates from the side of the battery, 
the following amounts of gassing were determined. Slight 
gassing was determined when occasional bubbles rose to the 
surface. This was about the same as when the battery was 
Standing. Light gassing occurred when there was a regular 
Stream of these bubbles rising to the surface. Moderate gas- 
Sing existed when considerably more activity was observed. 
This was accompanied on the charging curve by a noticeable 


increase in the slope. Heavy gassing was typified by a sheet 
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of small bubbles rising and a generally more violent condition 
at the surface of the electrolyte. This occurred when the 
slope of the charge curve reached its maximum and continued 
throughout overcharging when the voltage stopped increasing. 

To obtain equivalent charged conditions, the battery was 
overcharged when essentially no further increase in voltage 
occurred at constant current for about 15 to 30 minutes. Most 
of the energy into the battery at this point went into the 
hydrolysis of water, indicating that the charge reaction was 
essentially completed. The observation of the "coup de fouet" 
after this type of charge and its absense after a partial 
charge, adds some validity to the degree of completion of the 
reaction. 

An explanation of this phenomenon is that upon initiation 
of discharge, the electrolyte goes into a meta-stable state, 
super-saturated with PbSO,, until PbSO, crystallization sites 
Start to form on the plates. Once this occurs, the voltage 
recovers from its "dip" and proceeds in a smooth fashion. This 
recovery was observed to occur after about two minutes. The 
"coup de fouet" was not observed on a discharge following a 
partial charge, and was significant for discharges at higher 
discharge rates, 3.6 and 5.4 amps, which followed a full 
Charge. The magnitude of this dip was approximately O2.0L volt 

When the battery was new, Cell #2 had the lowest capacity. 


After several days of testing, while moving the experimental 
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set-up, Cell #3 was inadvertantly placed across a .12 resistor. 
After that occurred, Cell #3 had the lowest capacity and Cell 
#1 had the greatest. 

After all testing was completed, the following observa- 
tions of the battery were made. The positive plates and 
connectors had a dark brown color (artist's color "Burnt 
Umber"), which had a fairly rough appearance. There was some 
spalling on one edge of one of the plates in cell #3. 

The negative plates had a battle-ship grey color. One 
small piece of material was observed to have broken away from 
the plate. 

There was a fairly uniform coating of sediment in the 
bottom of the battery. This sludge consisted of two distinct 
materials. One was a light-colored material with a slightly 
pinkish tint. The other which uniformly covered the bottom 
of the battery, and was layered on top of the other sediment, 
was a dark brown color much the same as the color on the 
positive plates. The light-colored material may have been 
white rather than slightly pink because it had been observed 
by the author et al that the battery case itself seemed to cast 
the pink tint onto the interior components. It is believed 
that the light-colored material is actually predominantly 
PbSO,. 

Vinal (1955) states that overcharging produces sludging of 


a fine brown sediment, and that prolonged undercharging causes 





$7 
a fine white powder which is principally PbSO, to be deposited 
on the bottom. It is felt that the dark brown sediment is 
PbO2, which represents a loss of active material from the 
positive plates. The author notes that prolonged undercharging 


took place during life-cycle testing. 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 
The Shepherd (1965) equation 


Se Vglirg) - Voli) - Vi (q) (S21) 


V , (q) = activation polarization = A exp <i sa) 


(neglected due to the short time duration) 


Vo(i) = ohmic polarization = Rol (55) 


i (5.4) 





Vq(i-g) = diffusion polarization = K Oo - g 
proved to be an excellent starting point for establishing an 
electrochemical state model for lead-acid batteries. Using the 
lead-acid battery of this investigation, the author was able to 
get a fit to individual discharge curves which required only 
the modification of the ohmic polarization term (equation 5.3) 


to the form 


Vo(i,g) = (Rig + R,) i. (S75) 


in order to get an excellent fit. The results obtained from 
this fit suggested that a better fit for the family of curves 


Might be obtained by setting K*i constant. 
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The resulting modified term was 





Va(q) = Kee. (5.6) 


This was the form of the transport polarization suggested by 
Bode (1977). They also show that the use of the Peukert 


(1897) capacity model 
g(r) = er 7) (5.7) 


would predict the rapid drop of voltage at the end of discharge 
very well. 


The modification of the Va and Q terms seemed valid 


modeling for the family of discharges, but difficulties were 
apparent with the ohmic polarization term. There were no 
particular consistencies in the values of Ro and Ry for the 
different discharge rates. This may be due to the fact that 
the buildup of PbSO, at the plates occurs in a different manner 
for different discharge rates. Bode (1977) shows that for 
discharge rates of 6 mA/cm* and 1 mA/cm’ the distribution of 
PbSO, throughout the positive, capacity limiting plate, is 
uniform at full discharge. For current densities of 30 mA/cm* 


and 180 mA/cm’*, there is a lower PbSO, content at the center of 


the plate while the content near the surface is the same as 
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for the lower current densities. 

The cell under test had current densities of 3.4, 8.6, 
mie. 7, and 31 mA/em- for discharges of 0.6, 1.5, 3.6, and 5.4A 
respectively. This suggests that the same mechanisms were 
encountered in the battery discharge. Since this effect is 
clearly a function of current, a modification of Vo (i,q) other 
than that suggested in equation 5.5 may yield a better fit. It 
may have a form similar to Shepherd's original Va term, 
equation 5.4. 

Another consideration might be that since equation 5.4 was 
based on a small current density approximation, another look 
might be made toward applying equations 2.7 and 2.8 upon which 
the approximation was based. 

A very exciting area of future endeavor, and one which 
could prove very rewarding, would be the establishment of a 
model which would predict the family of discharge curves. 

Tremendous advancements in the lead-acid battery have been 
made since the time of Planté and this progress should continue 
steadily due to the importance of this device in myriad diverse 
applications. From the standpoint of the System Designer, the 
establishment of a good state model is a very exciting pros- 
pect. This would then allow for a determination of the optimum 
use of the battery, subject to the constraints of the applica- 
tion. 


Further effort in the direction of establishing a model 
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which would predict the family of discharge curves should 
include the effect of discharges at various rates during a 
Single cycle. Other effects such as pulsed discharge including 
both the frequency and duty cycle could then be added. 

The inclusion of the prediction of the change in internal 
characteristics during the secular life of the battery is 
virtually essential for an optimization. The objective in the 
design of almost any system is the optimization of a cost 
(possibly acquisition cost or life cycle cost) subject to 
various operational criteria. For a system such as an electric 
vehicle, a very significant cost is the repair cost of battery 
replacement. Therefore, it is necessary to be able to predict 
the lowest cost scheme for charge and discharge battery requir- 
ments and replacement frequency for various charge/ discharge 
schemes. To this end, further life cycle testing is recommend- 
ed. 

To produce applicable results, it is felt by the author 
that all cycles should be identical except those due to secular 
changes in the battery. A rather sophisticated control system 
needs to be developed. All charges and discharges should he 
Conducted at constant current. During charge, a significant 
period on overcharge is recommended to ensure that the state- 
of-charge is nearly identical prior to a discharge. Since 
the amount of overcharge cannot be determined by sensing the 


voltage level alone, some means needs to be devised to deter- 
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mine the slope of the voltage curve and to control the time 

on charge at this level. It is also recommended that a single 
cell be used during testing to eliminate any effects from the 
other cells. 

An attempt has been made throughout this thesis to draw 
conclusions on the suitability of the Shepherd (1965) equation 
as a state model for lead-acid batteries. Modifications to 
this equation were proposed, among them was the modeling of 
the ohmic polarization as a linear function of accumulated 
charge to obtain a good fit to a single discharge curve. These 
curves revealed the interesting result that the product KeI 
was nearly constant. Also, Peukert's equation was used to 
model the cell capacity. 

The difficulties encountered in the attempt to obtain an 
optimum state model have been evaluated and enumerated. How- 
ever, new discoveries about the lead-acid battery are being 
made, and new technology may supply the means to overcome the 
difficulties. Almost certainly, economic, ecological, socio- 
logical and technological problems with our current energy 
resources will dictate to society the necessity of assessing 


the place of the lead-acid battery in our future. 








93 


BIBLIOGRAPHY 


REFERENCE BOOKS 


Bode, H., Lead-Acid Batteries, John Wiley, New York, 1977. 


Kordesch, K.V., Batteries, Vol. 2, Marcel Dekker, New York, 


POT, 


Millman, J., and Taub, J., Pulse, Digital, and Switching 


Waveforms, McGraw-Hill, New York, 1965. 


ee 


Vinal, C.W., Storage Batteries, 4th ed., John Wiley, 


New York, 1955. 


REFERENCES 


Baikie, P.E., Gillibrand, M.I., and Peters, K., Electrochim, 


Acta, 17, 839 (1972). 


fase smi, | LeCthocectin. 4., 13; 857 (1966) . 


Bavytan, O.K., Bull. Acad. USSR (1946), 737. 








94 


eladstone, J-H. and Tribe, A., Nature, 25, 221, 461; 


Coie co 242,002; 27, 583, (1882-1883). 


Liebenow, C., Dissertation Gottingen 1905, 24. Elektrochem. 


4, 58-63, (1897). 


Pomkere, W., Elektr. Techn: 2.718 (1697), 287-288. 


Rabl, M., 4. Elektrochem. 42, (1936), 114-120. 


pemeoder, U., Electrotechn. 12 (1894), 587. 


Shepherd, C.M., J. Electrochem. Soc. 112 (1965), 657-664. 





75 
APPENDIX A 


NUMERIC DATA CORRESPONDING TO FIGURE 4-1. 
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Numeric data corresponding to Figures 4-4 and 4-5 
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Numeric data corresponding to Figures 4-6 and h-7. 
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Numeric Data corresponding to Figures 4-8 and 4-9. 


Coertficients for the model. 


meee table 4] 


0.6 Amp Discharge 1.5 Amp Discharge 





Ah Volt VOLE Dats Ah Vole Vole Blsese 
(Exp) (Eqn) (Exp) (Eqn) 
Ged foe 8 Za OS = Ole 1 Oo Oe4GS Pals? eyo aS Pe ed 
Peet 22060 2eS62 =C.C52 Sloe wacko fig SC eee © 
Se EC ek 0 Vee Seo Ce. 1 Le ies Padz2 cu ahe oe ence 
mee p tooo? (1.990 2.oonreg Lao SOO occ eee 
4,.8i bea? 9 bevee let he Zac 1699 See ee 
5240 leit 1.9639 Gea See Veg = Peer > Sed 
ee fC iene te ey Ce SFT Gs, Te Mette we. OG Salt 
GefD 12839 1-367 A By a7 fe oe ieee 30 cee 
Bele 1-810 ec. Oe &3 rae Vee be. Tee i Oe 
be 2S ees Ui 1efC9 ve Goo Ss tes ot Vet ta = Ogee 
6e3f 1.2/4 1¢/74& 0.008 te re, Pry oa gr a os i ee re Ce 
6236 iets leGs5 =~Oe 54 Geer Need. (pee ae Oem y a 
Bet tees t Peooo =—Oel 2s We ist Toe loo a = eS 
ae4 2? 164795 fe Ol) fs eG ely Ss 12660 View ee Se oe 
£244 eS © jen DU Cel Zt Bre ae. gm i a PA 
| Se eee. 8 Se 4.038. ei cs 





—_ >. 
Tea . fe 





a9 


Numeric data corresponding to figures 4-10 and 4-11 
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Numeric Data corresponding to Figures 4-14 and 4-15. 
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Numeric Data corresponding to figures 4-16 and 4-17 
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Numeric Data corresponding to Figures 4-19 and h-20. 
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Numeric data corresponding to Figures 4-21 and 4-22. 
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Numeric data corresponding to Figures 4-23 and -2h. 


COEFFICIENTS FOR THE MODEL 


SEE TABLE 4-3. 
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Numeric data corresponding to Figures 4-25 and h-26. 


COEFFICIENTS FOR THE MODEL 


SEE TABLE 4-3. 
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Numeric data corresponding to Figures 4-28 and 4-29. 
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Numeric Data corresponding to Figures 4-30 and 4-31. 


COEFFICIENTS FOR THE MODEL 
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Numeric Data corresponding to Figures 4-35 and 4-36 
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